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ANALYSIS OF THE LINEARIZED SUPERSONIC FLOW ABOUT POINTED 
BODIES OF REVOLUTION BY THE METHOD OF CHARACTERISTICS 

SUMMARY 


The method of characteristics applicable to the linear potential equation 
governing the supersonic flow over bodies of revolution at zero and small 
angles of attack is presented. For axisymmetric flow the equations used are 
those derived by Sauer and Heinz, whereas for flow due to angle of attack the 
equations are those derived by Oswatitsch and Erdmann. This analysis in- 
vestigates the limits of applicability of these methods. The flow field computa- 
tion was programmed in Fortran IV to be used for computing the axisymmetric 
and angle of attack flow fields about convex and concave parabolic ogives and 
cone-cylinder bodies of various body fineness ratios at low supersonic Mach 
numbers. The results of these computations were compared with J. L. Sims' 
exact method of characteristics, the linear theory of Reference 1, and, when- 
ever available, with experimental data. It was found that these methods are 
applicable to flow fields where the hypersonic similarity parameter, M m e, 
is less than about 0. 6. 


INTRODUCTION 

The linearization of the potential equation for supersonic flow about a 
body of revolution rests on the assumption that the disturbances produced in 
the flow by such a body are everywhere small. The solution of the potential 
equation may be obtained by either the method of singularities or the method 
of characteristics. The application of singularity methods to the present problem 
is well documented in the literature having originated with Von Karman and 
Moore [2] for axial flow, and Tsien [3] for bodies at angle of attack. The 
method of characteristics as applied to the linearized potential equation is not 
very well known, especially in this country. 

For flow fields where the potential equation can be linearized, there are 
several advantages in using linearized methods as opposed to exact methods. 

The speed of the linear methods and the simplicity of the equations involved are 
the primary advantages. Also, in the case of flow about cone-cylinder bodies, 
an important consideration is the ease with which the solution at the expansion 
corner may be obtained. In this situation, the exact method of characteristics 
involves a Prandtl-Meyer calculation in order to treat the flow field caused by 
the body slope discontinuity. In the linearized method of characteristics, the 



flow field calculation at the expansion corner is treated by the use of a "double" 
characteristic reflected from the corner. 

The linearized method of characteristics which is considered in this 
analysis is that of Sauer and Heinz [ 4] for axial flow, while for flow due to angle 
of attack the method is that of Oswatitsch and Erdmann [ 4] . The most important 
feature of these methods is that the particular combination of flow variables 
in the compatibility relations eliminates the singular behavior caused by the 
use of cylindrical coordinates so that these modified variables change quite 
slowly near the body. This in turn permits the use of a relatively large grid 
size which makes possible a very fast numerical calculation of the flow field. 

The nondimensionalization of the quantities used in the following treat- 
ment, where starred (* ) quantities are dimensional and unstarred quantities 
are dimensionless, is 



X* 

I* u* 

x L* 

u* 

u* 


V* 

R* (L * ) 

U* 

L* 

w* 


U* 

'-’oo 



The cylindrical coordinate system (x, r , ip) used in the following treat- 
ment is affixed with the body so that the x-axis is coincident with the axis of 
symmetry of the body (Fig. 1). 

TECHNICAL DISCUSSION 
The Nature of the Procedure 


We start with the linearized potential equation for supersonic flow 


«-l> 



a 2 # _ j_ a 2 $ 19$ 

9r 2 r 2 dip 2 r 9r 


( 1 ) 


2 



With G(x, r) as the perturbation potential of the axisymmetric part of 
the flow and F(x, r) as the perturbation potential of the flow due to angle of 
attack, we may write as a solution of equation (1) 


$ = cos a [x+ G(x, r)] + sin a cos ip F(x, r) . 


(2) 


The velocity components in terms of the potential <& are 

_ a$ _ a$ _ i a# 
u ax * v ar * w r dip ' 

Taking the inclination in the plane tp = 0, and introducing equation (2) 
into equation ( 1) , we get for the axisymmetric part of the flow 


,-- 2 « v 9fG 9fG i 8G _ n 
( oo " ) 0x 2 " 8r 2 " r 8r 


(3) 


and for flow due to angle of attack 


/ Twr2 _ U §!F_d?F_ldF±± 
00 9x 2 8r 2 r 9r r 2 


(4) 


are 


The axisymmetric disturbance velocity components divided by cos a 


u' = 


9G 

ax 


r' = 


8G 

8r 


and the angle-of-attack disturbance velocities divided by sin a cos tp are 


u" = 


8F 

ax ’ 


,M = 


8F 

8r 


To facilitate further derivation, the relations of equation (3) and (4) may be 
combined into a single relation 


... 8^ O <±> i OV , A 

( Mf, -1) ~ — T ~ _ 2 “ — 7 *■ ~T $ ~ 0, 

' °° ' a x z dr z r dr x l 


a z $ i a# 


(5) 


where 


$ = G for \ = 0 
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and 


$ = F for A = 1. 

The potential equation (5) can be written as 


9 

9u 

9v 

1 ,A „ 

r 

9x 

9r 

“ ~ v + 3" $ ~ 0 
r r 

9u 

9v_ 

0 


9r 

9x 



where 


( 6 ) 


/3 2 = M 2 -1, 

^ OO 


9# 

9x 


= u , 


9$ 

9r 


= v. 


This hyperbolic partial differential equation (for > 1 ) has two families 
of characteristic curves £ and 17 . The curves for £ = constant are referred to 
as left -running characteristic lines, while those for 77 = constant are called 
right-running characteristic lines. 

The following relations apply to these characteristic lines: 

| = x - /3r ( left-running lines) (7) 

17 = x + /3r (right-running lines) . 

From equations (1) through (7), the compatibility relations for axi- 
symmetric flow may be obtained. These relations. 


and 


9 ( rv* 1 

3 | 


9 ( rv* 1 
9tj 




( 8 a) 


( 8 b) 


are the fundamental relations on which the Sauer-Heinz linear characteristics 
method for axisym metric supersonic flow is based. They are used to obtain 
the values of u' and v’ along the characteristic lines described previously, 
which are just the Mach lines of the freestream. 
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The compatibility relations for flow due to angle of attack are obtained 
from equations (1) through (7) also. These relations 


and 



H 9(ru”) 
r d£ 



n 9(ru") 
r dr) 


(9a) 


(9b) 


are the relations on which the Erdmann-Oswatitsch linear characteristics 
method for flow due to angle of attack is based. They are used to determine 

F 

the values of ru" and v" + — along the characteristic lines, i. e. , the free- 
stream Mach lines. 


To solve the linearized potential equation (1) , it suffices to solve, 
in the case of axisymmetric flow, equations (8a) and (8b) for u' and v* and, 
in the case of flow due to angle of attack, equations (9a) and (9b) for u" and 
v". Knowing these values, any other desired flow property may be obtained 
at any point within the region of disturbed flow. 


Calculation of the Flow Properties in the Field 


To explain the calculation of the flow properties in the field, we will 
consider point 4 as an example (Fig. 1). The physical coordinates of the field 
points are obtained by solution of the left- and right-running characteristic lines 
simultaneously. Replacing equation (8a) by its finite difference form and ap- 
plying it along the right-running characteristic line 3, 4 we obtain 

(rv') 3 - (rv') 4 = - /3r^ (u 3 - u|). (10a) 

Similar application of equation (8b) along the left-running character- 
istic line 2, 4 yields 

( rv* ) 2 - (rv') 4 = /3r^ (u 3 - u 4 ) , (10b) 

where r^ is the average ordinate between the two points under consideration 
on the right- running characteristic line and r^ is the average ordinate between 
the two points under consideration on the left-running characteristic line. 
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Simultaneous solution of equations (10a) and (10b) yields for uj 


(rv*) 3 - (rv') 4 + /3(r^ u’ + r^ u^) 

«VV 


(ii) 


and for ( rv’ ) 4 


(rv') 4 



^ ~ U 3> 


( 12 ) 


Equations (11) and ( 12) give the values of the zero angle-of-attack flow 
properties, u' and rv', at field point 4. 


By the following identical procedure with equations (9a) and (9b), the 
relations which give the values of the angle-of-attack flow properties at point 
4 may be obtained. These relations are 


(ru") 4 = 


r. r 

_L_a 


IT v " + fj - ( v " + l ), ] + t [ r ; (ru ">l + V ru ">* ] 


+ V 


(13) 


and 



r „( V " + f l + ^( V " + f) g ~^[ (rU,,) 2- < rU,, ) 3 ] 

r + r 

t? I 


(14) 


Equations (11), (12), (13), and (14) give the values of the flow prop- 
erties at point 4 in terms of known values at points 2 and 3. The solutions of 
the other field points are obtained by a similar procedure. 


Calculation of the Flow Properties at the Body 

To explain the calculation of the flow properties at the body, we will 
consider point 5 (Fig. 1). The physical coordinates of the body point are 
obtained by simultaneous solution of the equation of the right-running charac- 
teristic line and equation of the body. 
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Consider a body contour given by r = h(x) (Fig. 1). Then the flow 

d$ 3$ 

direction © in the plane ip = 0, with velocity components u = — and v = — 
in accordance with equation (1), is given by the relation 


tan © = — 
dx 


9$ 

9r 

dx 


v 1 + v" tan a 
1 + u' + u M tan a 


( 15 ) 


Hence, there follows the exact boundary condition for axisymmetric flow: 


V 5 = (^ (1 + U 5 ) * (16) 

Equation (16) inserted into equation (15) gives, after brief manipulation, the 
exact boundary conditions for flow at angle of attack: 


v" = u" ~ • (17) 

Applying equation (8a) along 4,5, as in the previous section, we obtain 
(rv') 4 - (rv') 5 = - /3r n (uj - u^) . (18a) 

Simultaneous solution of equations (16) and (18a) yields 


0r u’+ (rv') 4 -r g , 

u 5 = 


*V^) S 


(18b) 


Equations (15) and (18b) give the values of the zero angle-of-attack 
flow properties at body point 5. 

Application of equation (9a) along 4, 5 yields 



(19) 
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Equation (19) shows that, in order to determine the value of v" at 
the body, an explicit knowledge of F along the body is necessary. The 
differential of F is 


dF(x, r) = — dx+ — dr = u" dx + v" dr. (20) 

If the step size is taken sufficiently small, the wall elements Ai (Fig. 1) can 
be considered straight lines, and we have 

Ax = cos © Ai and Ar = sin © Ai , (21) 

and from equation (17) , 

v" = u" tan ©. (22) 

As stated previously, the primary advantage of this method is that 
it permits the use of a relatively large grid size. The preceding restriction 
is not in contradiction with this fact since the method assumes that the rate 
of change in the slope of the body is small. 

Equations (21), which are difference equations, and equation (22) 
inserted into equation (20) yield 


AF (x, r) = 


u" Ai 2 

— Ai = u" i± — 

cos 0 Ax 


(23) 


Application of equation 23 along the body segment 2,5 gives F at the new body 
point 5 in terms of known properties at the previous body point 2 and u^' , 
which is still unknown: 


F = F + 77 [u"+ uTl . 

5 2 2 1 2 ~ J Ay 


Ai 2 
Ax 


(24) 


Combining equation (17) with equation (24) yields for (v" + ~) at point 5 


/ F\ /dh\ 

( v" + - ) = u" (— H- 
\ r / 5 \dx/ 5 


F 2 + « + »5> ^ 


(25) 


and equation (25) inserted into equation (19) yields, after some manipulation, 
u" at point 5 


8 


( 26 ) 



Equations (19) and (26) give the values of the angle -of-attack flow 
properties at body point 5. The solutions of all other body points are obtained 
in a manner similar to the one described above. 


The Method of Obtaining Starting Values 

It follows from the conditions for linearization of the potential equation 
that only pointed bodies of revolution can be treated by the procedure described 
in this analysis. 

Beginning at the apex of the body of revolution, the initial Mach line 
or left-running characteristic line is constructed. The inclinations of the 
characteristic lines are given by the relations 


Ar 1 

— = tan p = - 
Ax /3 


(left-running) 


and 


tan (-ju) = - j . (right-running) 

After choosing a suitable step size, the first point on the initial Mach 
line is located. In Figure i, this is point 1. The right-running characteris- 
tic line is drawn from this point and intersects the body at point 2. The portion 
of the body between the apex and point 2 is replaced by a cone formed by the 
chord line between these two points. Therefore, the flow in the triangular 
region formed by the apex, point 1, and point 2 is conical, thus enabling us to 
begin the characteristics solution. 

Because of assumed parallel flow at point i, we have u’ = u" = v* = 0, 

F 

F = r, and v" = 1. Therefore, we also have v" + — = 2. If we introduce the 

r 

easily proved geometrical relation 

_ h o / 1 \ 

% " 4 \ 3+ j3 tan / 
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into equation (18b) , we obtain for axisymmetric flow 


4 tan 2 1 9- 

U 2 ~ 1 + 3 /3 tan £ 0 + 4 tan 2 £■ Q 


( 27 ) 


Introducing the geometrical relation into equation (26) we obtain for flow due 
to angle of attack 


(1 + 3/9 tan S- ) sin 2 i? Q 

U 2 (1 + 3 /3 tan i? fl ) ( 1 + sin 2 t? Q ) + 4 /3 2 sin 2 t? 0 ’ 

where h 0 is the body radius at point 2 and 0 is the semi-apex angle. Equations 
(27) and (28) are independent of the body radius within their range of validity 

= ^0= CONSTANT 
dx 

and thus represent the approximate solution for conical flow. 

Cone-Cylinder Bodies 

In the solution of the flow field about cone-cylinder bodies by the methods 
described previously, we can see from the boundary conditions in equations (16) 
and (17) that a discontinuity in h, e. g. , at x = k in Figure 2, results in a dis- 
continuity in the calculation of the flow properties at that point. Therefore, it 
is necessary to calculate the flow properties at the body both upstream and 
downstream of the point of discontinuity, although the physical coordinates of 
the point are the same in each case. Because of the necessity of this double 
calculation at this point, the reflected characteristic from there produces a 
discontinuity in the calculation of the flow properties at each field point along the 
"double" characteristic. 

To explain the upstream and downstream calculation of the flow properties 
at the slope discontinuity on the body, we will consider points 5 and 6 of the 
configuration in Figure 2. The values of the flow properties upstream of the 
slope discontinuity are obtained by simpfy applying equations (16) , (18b) , (19) , 
and (26) at point 5. The zero angle-of-attack flow properties downstream of the 
singular point are found by reapplying equations (16) and (18b) at point 6; however, 
dli 

in this case, — , which is double-valued at the corner, assumes its downstream 
dx 

value. The relations which give the values of the angle-of-attack flow properties 
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downstream of the singular point cannot be used as they appear in equations 

A£ 2 

(19) and (26) due to the terms, — — and P 5 . It can be shown, however, that 

Ai 2 2 ^ X 

in the limit as Ax — 0, — - — 0 and that P does not change at the singular point; 

i. e. , F 5 = F e . Therefore, the relations which give the angle-of-attack flow 
properties at point 6 are with r g = r g = 


r 6 ( v " + f ) 5 + i 8< ru "> 5 - F S 



and 


ii _ it 

v 6 - u 6 



(30) 


The calculation of the field points on the left-running "double" character- 
istic from the corner follows the usual procedure outlined earlier, i. e. , the 
application of equations (11), (12), (13) and (14). The distinction of this cal- 
culation from any other lies only in the proper selection of base points on the 
left-running characteristic and in the case of the calculation of the downstream 
point (e. g. , point 10 in Fig. 2) and in the interpretation of r^. In this case 
is simply r 9 = r 10 since the points 9 and 10 are physically coincident. 

It is desirable in the cone-cylinder case to have several intermediate 
points on the right-running characteristic which intersects the expansion corner 
so that the flow field downstream of the corner will be as accurate as possible. 
For this reason we choose the initial step on the cone, i. e. , the two point 
conical approximation as given by equations (27) and (28) , to be a fraction of 
the total cone length. 


Pressure Relations and Force Coefficients 

According to Oswatitsch and Erdmann [4] , the pressure coefficient at 
angle-of-attack a is given by the relation 

(1!^ = — 2 (u' + u" a cos ip) + ex 2 

- v' 2 + /3 2 (u 1 + u" a cos ip) 2 


(l-^f sin 2 ip) 


(31) 
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For small angles of attack, the second order terms in a may be neglected, and 
we obtain 

C = C + C a cos ip, (32) 

P Pi P 2 

where C = /3 2 u' 2 - 2u' - v' 2 is the pressure coefficient at zero angle of 
Pl 

attack and C = 2u" (/3 2 u" - 1) . 

P2 

The normal force coefficient is given by 


sr- 


2 

S(L) 


L 7T 

I f 


o o 


Cp cos ip h d ip dx . 


(33) 


Introducing equation (32) into equation (33) and differentiating with respect 
to a, we get the normal force coefficient slope 


C 


N 


a 


d ( c N ) 


TT 

S(L) 


L 

f C • h • dx . 

i P2 


(34) 


It follows immediately from equation (34) that the local normal force variation 
is given by 


d ' C N > 
a 

dx 


TT 

S(L) 


C • h. 
P2 


RESULTS AND DISCUSSION 


The Sauer-Heinz and Oswatitsch-Erdmann methods were programmed in 
Fortran IV to solve flow fields varying in body shape, body fineness ratio, and 
Mach number. To investigate the limits of applicability of these methods, the 
results of these computations were compared with J. L. Sims' exact method of 
characteristics* , the linear theory of Reference 1, and experimental data. 


Sims, J. L. : Results of Method of Characteristics Analysis at Small 
Angle of Attack (Unpublished). Aero-Dynamics Division of NASA, 
Marshall Space Flight Center, 1965 
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The approach taken by Sims is based upon a small angle-of-attack pertur- 
bation superimposed on a nonlinear axisymm§tric flow field which is determined 
by the method of characteristics (eqs. 26-15, 26-16, 26-17, and 26-18 in Ref. 5). 


For comparison of the zero angle-of-attack flow properties computed 
by the Sauer-Heinz method, the pressure coefficient was computed for convex 
and concave parabolic ogives. The equations of these ogives are 

h(x) = 7 (x + x 2 ) concave ogive 

A 


and 


h(x) = e (2x - x 2 ) . concave ogive 

The values of the body fineness ratio, e, used in these computations were 

0. 05, 0. 10, 0. 15, and 0. 20, and the Mach numbers were 1. 5 and 3. 0. These 
results are given in Figures 3 through 18. The pressure coefficients were 
also computed for 5-degree and 10-degree cone-cylinder bodies at Mach . 
numbers of 1. 5 and 3. 0 and an ogival body whose configuration is given by the 

equation h(x) = irx - 77 x 2 [1} at Mach numbers of 1. 84, 2. 45, and 3. 01. 

3 36 

These results are given in Figures 19 through 25. 

For comparison of the angle-of-attack flow properties computed by the 
Oswatitsch-Erdmann method, the normal force coefficient slope was computed 
and plotted as a function of Mach number for the convex and concave ogives and 
the 5-degree and 10-degree cone-cylinders. These results appear in Figures 
26 through 31. Also, the local normal force variation was computed for the con- 
vex and concave ogives at Mach numbers of 1. 5 and 3. 0 and the ogival body [1 ] 
at Mach numbers of 1. 84, 2. 45 and 3. 01. These results appear in Figures 32 
•through 42. 

For comparison of the Oswatitsch-Erdmann method with experimental 
data, the local normal force variation was computed for a convex parabolic 
ogive with body fineness ratios of 0. 066, 0. 1, and 0. 1667 at Mach numbers of 

1. 59, 3. 0, and 4. 25 and also for a 10-degree cone-cylinder at Mach numbers of 
1. 99 and 2. 44. These results appear in Figures 43 through 49. 

Figures 3 through 18 show that the Sauer-Heinz method agrees closely 
with Sims' method for body fineness ratios of 0. 15 or less; however, for a 
body fineness ratio of 0. 20 they begin to deviate considerably, especially at the 
higher Mach number. These two methods compare very well for the 5-degree 
and 10-degree cone-cylinders also. In view of the assumptions made in the 
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linearization of the potential equation, this is exactly the trend that one would 
expect. This same trend can be observed when comparing the Oswatitsch- 
Erdmann method with Sims' method in Figures 26 through 39. Although the 
range of values for the body fineness ratios and Mach numbers in the experi- 
mental data is somewhat limited, the trend seems to be present there also. 
Figures 19 through 21 and 40 through 42 show that the comparison between these 
methods and the linear theory of Reference 1 is very good except at the apex of 
the body. This disagreement there could be attributed to the size of the initial 
step-length used in the conical approximation flow region, from which the start- 
ing values are obtained. 


CONCLUDING REMARKS 


The truncation error in the numerical solution was studied by system- 
atically varying the field mesh width. A constant step-length in the entire 
field was used. Step-lengths in the range 0. 0025 -*■ 0. 0200 were tried to 

increments of 0. 0025. The truncation error was found to be sufficiently small 
with a step-length of 0. 01, which was then used to compute all of the solutions 
for this report. Using this step-length, the average time for a complete run 
for both zero and nonzero angle-of-attack solutions at one Mach number was 
about thirteen seconds on the IBM 7094 digital computer, 

A careful study of the comparisons of the Sauer-Heinz and Oswatitseh- 
Erdmann linear methods given in the figures shows that these methods are 
applicable to flow fields where the hypersonic similarity parameter, Mm e, 
is less than about 0. 6 


George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 
Huntsville, Alabama, April 19, 1966 
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FIGURE 4. PRESSURE DISTRIBUTION ALONG A CONVEX PARABOLIC OGIVE WITH 
BODY FINENESS RATIO OF 0. 05 AT = 3.0 
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FIGURE 7. 


FIGURE 8. 
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PRESSURE DISTRIBUTION ALONG A CONVEX PARABOLIC OGIVE WITH 
BODY FINENESS RATIO OF 0. 15 AT M*, = 1. 5 
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Sims' Exact Method of Characteristics 
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PRESSURE DISTRIBUTION ALONG A CONVEX PARABOLIC OGIVE WITH 
BODY FINENESS RATIO OF 0. 15 AT M*, = 3. 0 
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Sims' Exact Method of Characteristics 
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FIGURE 11. PRESSURE DISTRIBUTION ALONG A CONCAVE PARABOLIC OGIVE WITH 

BODY FINENESS RATIO OF 0.05 AT M =1.5 
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FIGURE 12. PRESSURE DISTRIBUTION ALONG A CONCAVE PARABOLIC OGIVE WITH 

BODY FINENESS RATIO OF 0. 05 AT M„ = 3. 0 
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FIGURE 13. PRESSURE DISTRIBUTION ALONG A CONCAVE PARABOLIC OGIVE WITH 

BODY FINENESS RATIO OF 0. 10 AT M = 1. 5 


FIGURE 14. 









FIGURE 15. 


FIGURE 16. 
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PRESSURE DISTRIBUTION ALONG A CONCAVE PARABOLIC OGIVE WITH 
BODY FINENESS RATIO OF 0. 15 AT = 1.5 



PRESSURE DISTRIBUTION ALONG A CONCAVE PARABOLIC OGIVE WITH 
BODY FINENESS RATIO OF 0. 15 AT M m = 3. 0 
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FIGURE 17. 


FIGURE 18. 
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PRESSURE DISTRIBUTION ALONG A CONCAVE PARABOLIC OGIVE WITH 
BODY FINENESS RATIO OF 0. 20 AT M*, = 1.5 
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PRESSURE DISTRIBUTION ALONG A CONCAVE PARABOLIC OGIVE WITH 
BODY FINENESS RATIO OF 0. 20 AT = 3. 0 
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— Linearized Method of Characteristics 
O Linear Theory, Reference i 
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PRESSURE DISTRIBUTION ALONG A BODY CONFIGURATION DESCRIBED 
BY THE EQUATION, h(x) = x/3 - 10x 2 /36, = 1. 84 
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FIGURE 20. PRESSURE DISTRIBUTION ALONG A BODY CONFIGURATION DESCRIBED 
BY THE EQUATION, h(x) = x/3 - 10x 2 /36, = 2. 45 
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PRESSURE DISTRIBUTION ALONG A 5 DEGREE CONE PLUS 
CYLINDER BODY, = 3. 0, e » . 05 
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FIGURE 24. PRESSURE DISTRIBUTION ALONG A 10 DEGREE CONE PLUS 

CYLINDER BODY, = 1. 5, e = . 10 
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FIGURE 27. EFFECT OF MACH NUMBER AND THICKNESS RATIO ON NORMAL 
FORCE COEFFICIENT SLOPE FOR A CONVEX PARABOLIC OGIVE 



FIGURE 28. EFFECT OF MACH NUMBER AND THICKNESS RATIO ON NORMAL 
FORCE COEFFICIENT SLOPE FOR A CONCAVE PARABOLIC OGIVE 




Linearized Method of Characteristics 

• Sims 1 Exact Method of Characteristics 



Mach Number 


FIGURE 31. EFFECT OF MACH NUMBER ON THE NORMAL FORCE COEFFICIENT 
SLOPE FOR A 10 DEGREE CONE PLUS CYLINDER BODY 



FIGURE 32. LOCAL NORMAL FORCE VARIATION FOR CONVEX AND CONCAVE 
OGIVES OF FINENESS RATIO 0.05 AT A MACH NUMBER OF 1.5 








FIGURE 35. LOCAL NORMAL FORCE VARIATION FOR CONVEX AND CONCAVE 
OGIVES OF FINENESS RATIO 0. 10 AT A MACH NUMBER OF 3. 0 



FIGURE 36. LOCAL NORMAL FORCE VARIATION FOR CONVEX AND CONCAVE 
OGIVES OF FINENESS RATIO 0. 15 AT A MACH NUMBER OF 1. 5' 
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FIGURE 37. LOCAL NORMAL FORCE VARIATION FOR CONVEX AND CONCAVE 
OGIVES OF FINENESS RATIO 0. 15 AT A MACH NUMBER OF 3. 0 



FIGURE 38. LOCAL NORMAL FORCE VARIATION FOR CONVEX AND CONCAVE 
OGIVES OF FINENESS RATIO 0. 20 AT A MACH NUMBER OF 1. 5 
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FIGURE 39. LOCAL NORMAL FORCE VARIATION FOR CONVEX AND CONCAVE 
OGIVES OF FINENESS RATIO 0. 20 AT A MACH NUMBER OF 3. 0 
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FIGURE 40. LOCAL NORMAL FORCE VARIATION ALONG A BODY CONFIGURATION 
DESCRIBED BY THE EQUATION, h(x) = x/3 - 10x z /36, = 1. 84 
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FIGURE 43. COMPARISON BETWEEN THEORY AND EXPERIMENT OF LOCAL NORMAL FORCE 
VARIATION FOR A CONVEX OGIVE OF FINENESS RATIO 0. 066 AT A MACH NUMBER OF 1. 59 

Linearized Method of Characteristics 

Sims' Exact Method of Characteristics 


©Experiment, Reference 7 



X 


FIGURE 44. COMPARISON BETWEEN THEORY AND EXPERIMENT OF LOCAL NORMAL FORCE 
VARIATION FOR A CONVEX OGIVE OF FINENESS RATIO 0. 1 AT A MACH NUMBER OF 3. 0 
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FIGURE 45. COMPARISON BETWEEN THEORY AND EXPERIMENT OF LOCAL NORMAL FORCE 
VARIATION FOR A CONVEX OGIVE OF FINENESS RATIO 0. 1 AT A MACH NUMBER OF 4. 25 
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FIGURE 46. COMPARISON BETWEEN THEORY AND EXPERIMENT OF LOCAL NORMAL FORCE 
VARIATION FOR A CONVEX OGIVE OF FINENESS RATIO 0. 1667 AT A MACH NUMBER OF 3. 0 




FIGURE 47. COMPARISON BETWEEN THEORY AND EXPERIMENT OF LOCAL NORMAL FORCE 
VARIATION FOR A CONVEX OGIVE OF FINENESS RATIO 0. 1667 AT A MACH NUMBER OF 4. 25 



FIGURE 48. LOCAL NORMAL FORCE VARIATION ALONG A 10 DEGREE CONE PLUS 

CYLINDER BODY, M^ = 1. 99, e = . 10 
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FIGURE 49. LOCAL NORMAL FORCE VARIATION ALONG A 10 DEGREE CONE PLUS 

CYLINDER BODY, M* = 2.44, e = . 10 
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"The aeronautical and space activities of the United States shall be 
co?iducted so as to contribute ... to the expansion of human knowl- 
edge of phettomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof ” 

— National Aeronautics and Space Act of 1958 
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